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^ '. Abstract 

•/^ I We show that the entanglement of an initially frequency-entangled photon pair propagating 

^H ■ through a disordered scattering medium can be read off in transmission from the coincidence 

a ■ 

I . counting rate. The very same quantity also encodes the Thouless time of the medium. 



C^ 



o 



> 

oo 

o 

(N 

O 

o 



X 



PACS numbers: 42.25.Dd, 42.50.Lc, 42.65.Lm 



I. INTRODUCTION 

The study of speckle patterns, irregular interference structures produced by light propa- 
gating through an inhomogeneous medium, is of fundamental importance in the understand- 
ing of coherent wave propagation in disordered media [1]. While, at first sight, a speckle 
pattern looks perfectly random, closer analysis reveals a more complex picture. In the 
backscattering direction, for instance, the light intensity is enhanced, due to the construc- 
tive interference of pairs of counter-propagating multiple scattering paths. This is known 
as the coherent backscattering effect [2]. In transmission, speckle patterns exhibit long- and 
infinite-range correlations between distant spots. The latter are, in particular, at the origin 
of universal conductance fiuctuations in electron transmission across mesoscopic conductors 
[3] . Once again, these highly nontrivial correlations originate from the interference between 
pairs of scattering paths which persists under disorder average, and have been studied ex- 
tensively during the last decades [1, 4]. 

While our knowledge of the statistics of classical waves in disordered media is today fairly 
well developed, little is known about the statistical properties of nonclassical light propa- 
gating in such media. Recently, an important step forward was made with the investigation 
of the interplay between quantum fiuctuations of light, on the one hand, and the classical 
fiuctuations induced by a scattering medium with static disorder [5, 6] or moving scatterers 
[7] , on the other hand. In this scenario, a given quantum state of light enters a disordered 
medium, and signatures of the quantum statistics are found in the photon number variance 
analyzed upon transmission. Later, the probability distribution of the coincidence counting 
rate, which defines a "two-photon speckle", was analyzed on output of a disordered medium 
[8]. Two-photon speckles are relevant for studying entangled states of light, and were re- 
cently observed experimentally [9], but interferences between different scattering paths were 
assumed to vanish under disorder average. 

Interference phenomena with nonclassical light in disordered media have so far only been 
addressed in the framework of one-dimensional tight-binding models [10] or in the transport 
of quantum states with independent photons [11]. In our present contribution, we consider 
the propagation of a pair of spectrally entangled photons in an open, disordered medium 
with three-dimensional disorder. We show that the emerging interference effects bear all 
the information on the spectral entanglement between the two photons, and that, more 



surprisingly, information on the dynamical properties of the disordered medium is also in- 
scribed in the interference signal. Finally, we suggest a method to access this information 
experimentally. To do so, we analyze the photon coincidence counting rate of a photon pair 
transmitted through a disordered waveguide, within a continuous-mode approach (Sees. II, 
III). We show that interference contributions to the coincidence rate depend on two time 
scales, a time t^ characteristic of the spectral, and, in particular, the entanglement properties 
of the photon pair, and the Thouless time tD, i.e. the typical time it takes for a photon to 
traverse the disordered sample by diffusion [12]. In Sec. IV, we examine how the interference 
contributions depend on these characteristic times. Finally, in Sec. V we propose a way 
to measure these two times experimentally, by introducing a small delay between the two 
photons before they enter the medium. In this situation, we show that interference leads to 
a peak in the coincidence counting rate, whose size gives direct access to t,, and to. 

II. MODEL 



Our model is depicted in Fig. 1. A source produces pairs of photons which are scattered 
from a disordered waveguide of length L and cross-section A. The mean free path, i.e., the 
average distance between two consecutive scattering processes, is assumed to be small, i ^ 
L, such that each photon experiences multiple scattering before being collected on output by 
two detectors. In the most general situation considered in this paper (Sec. V), one of the two 
photons is delayed by a time 6t with respect to the other. If we assume perfectly reflecting 
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FIG. 1: (color online). Two photons, incident in spatial modes a and a', are scattered from a 
disordered waveguide of length L, with a mean free path i <^ L. The coincidence counting rate 
between two outgoing modes b and b' is analyzed by photodetectors. The input-output relation (2) 
takes into account the coupling to all unoccupied (vacuum) modes a and /3, which are depicted by 
dashed arrows. 



boundary conditions at the outer surface of the waveguide, only a finite number 2N{u) = 
4 X k{ujYA/{4TT) of transverse spatial modes are supported, where k{u) is the wave vector 
at frequency u. Note that N{u) accounts for two possible propagation directions and two 
polarization states of the electromagnetic field. Another important quantity characteristic 
of the disordered medium is the dimensionless conductance g = {A/3)Ni/L. g is the single 
parameter that controls the transport properties across the disordered waveguide, from the 
regime of Anderson localization {g ^ I) to the regime of diffusive multiple scattering g > 1 
[13]. We will focus on the case g > 1 hereafter. 

An incoming photon pair is described by a two-photon wave-packet, within a continuous- 
mode formalism, with the electromagnetic field quantized in the usual way [14]. The two- 
photon state is written in the general form {tp) = J dudu'S{u,u')al{uj)al,{uj')\vac), where 
the continuous-mode operator dl^{uj) creates a photon with frequency u in the spatial mode 
a. The two-photon wave-packet amplitude S{u,u') contains all spectral information about 
the photon pair. In analogy with [15], we choose 

1/2 



S{u,u') 



1 



2n,/T^/tl 



exp 



-{id - LUo) " (^' " ^o) + 2r7(cj - Wo)(^' - Wo) 






4^1 - nVtl 

(1) 

which provides a simple and general model of a two-photon wave-packet, characterized by a 
central photon frequency Wqj and a single-photon coherence time tc- Experimentally, such a 
two-photon wave-packet can be produced by means of parametric down conversion, in which 
two entangled photons originate from a nonlinear crystal irradiated by a "pump" photon. 
The dimensionless parameter r] G ] — 1, 1[ determines the frequency entanglement of the 
two down-converted photons, and depends on the coherence time of the pump photon, the 
length of the crystal, and the phase velocities of the pump and the down-converted photons 
in the crystal [16]. Different types of entanglement are distinguished by the cases r^ < 0, 
r] > and r] = 0. These correspond to anticorrelated, correlated and uncorrelated photons, 
respectively. The amplitude |?7| gives the "strength" of entanglement, the limits r^ — ?> ±1 
corresponding to fully entangled states. Fully anticorrelated photons (r^ — )■ —1) correspond 
to SibJ^bj') oc bibj + w' — IbjQ) and are produced by standard parametric-down conversion 
with a monochromatic pump [17]. The possibility to generate fully correlated [j] — )■ 1), 
for which S'(a;,w') oc (5(w — w'), or completely uncorrelated [ri = 0) photons, for which S 
can be written as the product of two independent wave-packets, has also been discussed 



[16, 18]. Recently, an original experimental method was devised to produce photon pairs 
with arbitrary values of f], with the help of a superlattice of nonlinear crystals [19] . Let us 
finally note that the Gaussian form of Eq. (1) is chosen for convenience, since the results 
derived herafter do not qualitatively change for slightly different shapes of S. 

In order to study the transport properties of the photon pair in the disordered medium, 
we make use of standard input-output relations between incoming and outgoing annihilation 
operators [20]: 

dbiuj) =^tabiuj)aai(^) + ^ rf3b{uj)ap{uj), (2) 

Q=l l3=N{uj)+l 

where tab and r/35 are transmission and reflection coefficients from the incoming modes a and 
/3, respectively, to the outgoing mode b. Note that Eq. (2) fully accounts for the coupling to 
all vacuum states a ^ a,a' and P ^ b,b' (see Fig. 1). Creation and annihilation operators 
obey the usual bosonic commutation relations [ai(cj), a](w')] = 5ij5{uj — cu'). 

III. COINCIDENCE COUNTING RATE FOR St = 

In this section, we assume that there is no time delay between the incident photons. We 



first examine the mean photon number (lb) = L dt {dl{t)db{t)) registered by a photodetector 
in a given outgoing mode b, during a sampling time T. This quantity gives the probability 
to obtain a "click" from the photodetector, meaning that it collected one photon of the pair. 
(■ ■ ■ ) refers to quantum mechanical averaging over the two-photon state, and ^^^ to classical 
ensemble averaging. These two types of averaging highlight that the photon number fluctu- 
ates due to both, the quantum nature of the light, and the stochastic properties of the disor- 
dered medium. The time-dependent operators db{t) = {1/\/2tt) J^"_^^ L dcudb{uj)exp{—iujt) 



(JQ— Aa;/2 

are defined via a Fourier integral over the bandwidth Au ^ uq of the photodetector around 
the central photon frequency Uq. Au is typically much larger than the spectral bandwidth 
1/tc of individual photon wave-packets. In addition, since usual photodetectors are not able 
to resolve photon wave-packets, the detection time T exceeds any other time scale of the 
problem, such that we have 

1/T < 1/tc < Aw < uo. (3) 

This separation of time scales allows to extend the range of frequency and time integrations 
from —00 to 00, and to approximate the number of transverse spatial modes by N{u) ~ 



A^(a;o). We further assume that the properties of the disordered medium do not change at 
the scale of Au, such that quantities hke the mean free path are constant upon integration. 
Under these premises, the input-output relations (2) and standard diffusion theory [4] lead 



to (Jfe) = {A/N)(i/L). Note that (/{,) is independent of t^ and rj, i.e. the mean photon 
number does not contain any spectral information on the two-photon state [8]. We therefore 



turn to the coincidence number {hh') between two outgoing modes b and b' (with possibly 
b = b'), which will be the figure of merit from now on. This quantity is defined as 

{hh') = dt dt'{dl{t)al{t')dbit')db{t)). (4) 

Jo Jo 

After ensemble averaging, this turns into 
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with Eqs. (2) and (3). The kernel Kaa'hh' accounts for fourth-order interference between 
photon path amplitudes multiply scattered by the disorder. These interferences are described 
by combinations of products of four transmission coefficients, averaged over disorder, and 
represented by Feynmann diagrams [4]. In the diffusive limit g > 1 considered here, Kaa'bw 
can be expanded in powers of 1/ g. Performing this expansion up to the first order in 1/ g, 
the coincidence counting rate R reads 



(h) X ih') 2 1 3g \tnj [ \tDj Sg \tnj 



(6) 



{h) X {h') 
which is the main result of the present contribution. In Eq. (6) we introduce two char 



acteristic times, t^ = t^ \/\ — rj and the Thouless time tn = L /ij D), with D being the 
photon diffusion coefficient. Whereas t^ contains the spectral information on the photon 
pair (through the coherence time t^. and the entanglement parameter 77), t£, is a quantity 
characteristic of the disordered medium, giving the typical time it takes for a photon to 
diffuse through the sample. The Thouless time was initially introduced by Thouless in the 
context of electronic conduction in metals [12, 21], and is, in essence, a quantity character- 
izing the dynamics of the disordered medium, due to its proportionality to the inverse of 
the diffusion coefficient [22-24]. The dependence of R on t^ and t£, is discussed in detail in 
the next section. Here we first comment on the physical content of Eq. (6): its right-hand 
side contains four terms. The first, constant term describes the independent propagation of 
the two photons in the disordered medium, and corresponds to the diagram in Fig. 2a. The 
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term denoted C^^^ is depicted in Fig. 2b. It results from a fourth-order interference, said 
to be local, because the exchange of photon path amphtudes (the "crossing" ) occurs close 
to the boundary of the medium. The two terms oc C*^^^ and oc Sbb'C^^^ in Eq. (6) encode 
another type of fourth-order interference, in which photon path amplitudes change partners 
at some point in the multiple scattering sequence. They are represented by the diagrams 
of Fig. 2c and 2d, respectively [25]. Note that here the photon path amplitudes recombine 
and then propagate again after the crossing, such that the underlying interference process 
is nonlocal. 
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FIG. 2: Leading order diagrams that contribute to the kernel Kaa'bb', and give rise to the four 
terms in the right-hand side of Eq. (6). In each diagram, the two parallel lines connected by 



dotted "ladders" symbolize averages of products of transmission coefficients, tij{uj)ti'j'{u}'), with 
incoming i,i' = a or a' and outgoing j,j' = b or b' directions [25]. 



We now discuss the structure of Eq. (6). When the coincidence counting rate is measured 
in one outgoing mode, b = b', such that R ~ l/2[l+C^^^ + {2/3g)C^^'^] ~ 1/2[1+C^% and the 
local interference is the main interference contribution to R. On the contrary, when b ^ b', 
R :ii 1/2 [1 + {2/3g)C^'^^] and the nonlocal interference is the main interference contribution 
to R. Therefore, C^^' is a "short-range" interference effect unlike C^^\ which is spatially 
long-range. These properties are reminiscent of classical optics in disordered media, where 
local interferences explain the granularity of speckle patterns, whereas nonlocal interferences 



are responsible for long-range correlations in the speckle pattern [1]. 

IV. INTERFERENCE, SPECTRAL INFORMATION AND THOULESS TIME 

The analysis of the previous section reveals that interference, i.e., terms proportional to 
C^^^ and C*-^-* in Eq. (6), depend on the ratio trj/tjj. Therefore, they contain information 



about the spectral properties of the pair, through the characteristic time t^ = td ^/T^~rj, 
and information about the dynamics of the disordered medium, through the Thouless time 
t£). In this section, we analyze more closely the dependence of C*^^-* and C-^^ on this ratio. 
As shown in Fig. 3, both C^^'> and C^^^ increase with tr,/tD, before they saturate at the 




FIG. 3: (color online). Interference contributions C^^' and C'^^ to the coincidence counting rate R 
[see Eq. (6)], plotted as a function of tr^/to- Both saturate at unity for tri/to ^ 1- 

asymptotic value C-^^ ^ C-^^ ~ 1, for t^ ^ to. In this limit, the disordered medium is 
insensitive to the spectral properties of the photon pair, and the coincidence rate reaches 
a constant value. This is the monochromatic limit, in which propagation of the pair can 
be treated within a single-mode approach [5, 11]. According to the definition of t.^, the 
increase of C-^^ and C*-^-* can be induced in two different manners: by an increase of the 
coherence time tc of individual photons, or of the entanglement parameter rj. While the first 
case simply means that interference between multiply scattered photons can only occur if 
the individual photon wave-packets overlap {i.e., if tc is not too small), the second is less 
obvious. To clarify the role of entanglement in the interference terms in Eq. (6), we show 
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FIG. 4: (color online). Functions C^^> (solid curves) and C*-^-* (dashed curves) plotted as a function 
of the entanglement parameter r/, for tc/^D = 0.1 (red curves), tc/^D = 1 (green curves) and 
^d^D = 8 (blue curves). All curves saturate at unity for strongly correlated photons (r? — )• 1). 

in Fig. 4 C^^' and C^^-* as a function of the parameter ?], for three different values of tc/tr). 
All curves increase with rj and saturate at unity when rj -^ 1. This can be understood as 
follows: when 77 < 0, a given frequency component of one photon tends to be entangled with 
a different frequency component of the other (quantum anticorrelation) , and therefore does 
not interfere with it, leading to a reduction of the interference signal after averaging over 
different realizations of the disorder. Conversely, when 77 > 0, a given frequency component 
of the first photon tends to be entangled with the corresponding same frequency component 
of the other (quantum correlation) , and therefore interferes constructively with it, leading 
to an enhancement of the interference signal after averaging. From Fig. 4, we note, however, 
that the behavior of anticorrelated and correlated photons is not symmetric. In particular, 
in the limit r] — )■ 1 of perfectly correlated photons, interference is fully constructive, and 
photons behave like monochromatic ones, even if the coherence time tc is small [15]. 

V. COINCIDENCE COUNTING RATE WITH DELAYED PHOTONS 



We have shown in the previous section that interference contributions to the coincidence 
counting rate (6) carry information (encoded in the ratio t^Jtj:,) on the entanglement prop- 
erties of the photon pair and on the dynamics of the disordered medium. We now propose 



an experimental scheme giving individual access to these two characteristic times. For this 
purpose, we assume that one of the photons arrives with a small delay 5t with respect to 
the other (see Fig. 1). For 5t ^ 0, the coincidence counting rate (6) becomes 
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(7) 

i.e., all interference contributions acquire a dependence on 6t. This dependence originates 
from the fact that the occurence of interference is only possible if the time intervals during 
which the two photons propagate in the medium overlap. C-^^ and C-^^ are shown in the 
main plot of Fig. 5, as a function of 5r/t^, for several values of tr^/tn- We recall that C*^^^ 
is the dominant interference contribution to R when b = b' (since i? ~ 1/2[1 + C-^^]), and 
C^"^^ when b ^ b' (since R ^ 1/2[1 + {2/3g)C^'^^]). Two different regimes can be identified 

C(i/2)(5r/t^) ^ 




2 2 

time delay (5r/t^ 



FIG. 5: (color online). Main plot: functions C'^-* (solid curves) and C^^-* (dashed curves) plotted 
as a function of the time delay Sr/tri, for trj/to = 0.4 (red curves), ty^/to = 2 (green curves) 
and t^/tD = 7 (blue curves). Inset: mean square size (a/t^)^ of the functions C^^>{6T/t^) and 
C^'^'{5T/t^), as a function of t^/tj:) (cj^ is defined in the text). For t^/tj:) ^ 1, cr ~ t^, whereas for 
tri/to ^ 1, o- ~ to- 

from Fig. 5: first, when t^ ^ to, C^^^ ~ C-^^ ~ exp[— ((5r/t^)^] are peaked functions of 
6t, and the mean square size a^ = J^^6T^C^^/^\6T)d{6T)/ J^^C^^/^\6T)d{6T) ~ 1^/2, as 
obvious from the asymptotic behavior in the inset of Fig. 5. In this limit, measurement of 
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R for different time delays makes the joint coherence time t^ of the photon pair accessible. 
In the opposite limit t^, ^ tu-, C*-^-* and C*-^^ are broad functions of 5r, with a reduced 
visibility. The inset of Fig. 5 clearly shows that o" oc t/) for both C^^\5t) and C'^'^\5t) in 
this limit. The coincidence counting rate then gives access to the Thouless time, and thus 
to the diffusion coefficient. 

In the diffusive regime considered here, nonlocal interference [which accounts for the C'-^-* 
terms in the coincidence counting rate (6)] is of weak magnitude, and is, therefore, more 
difficult to measure than local interference. When h ^ h' however, R ^ 1/2[1 + {2/'ig)C^'^^\ 
[see Eq. (7)], and the observation of nonlocal interference may be greatly facilitated by 
taking advantage of its dependence on 5r, which allows us to distinguish it from the constant 
background of 1/2. An interesting problem concerns the regime where gf ~ 1, i.e., when 
Anderson localization sets in. It was shown recently that nonlocal interferences persist in 
this regime [1 1] , where they are of the order of unity and acquire a dependence on the ratio 
of the sample size L and on the localization length. However, in the context of the present 
work, their dependence on the spectral properties of the photon pair and on the dynamical 
properties of the medium is less clear, in particular since the Thouless time is no longer a 
relevant quantity when (? ~ 1. 

VI. SUMMARY AND CONCLUDING REMARKS 

We have calculated the coincidence counting rate that results from the propagation of 
a pair of entangled photons through a disordered medium, and we have analyzed its de- 
pendence on the fourth-order interference that occurs between multiple scattering photon 
trajectories. We have shown that these interferences convey information on the spectral en- 
tanglement of the pair and on the Thouless time of the disordered medium. Finally, for the 
purpose of measuring this information, we have proposed an experimental scheme in which 
the photons of the pair are incident on the disordered medium with a time delay. In classical 
optics, the measurement of the Thouless time usually requires a dynamical experiment, such 
as time-resolved transmission of short pulses [22, 23] or spectral correlation mesurements 
[24]. In the method we propose, no dynamics are involved, only the counting of photon 
coincidences in the framework of a quantum optics experiment. 

Let us conclude with a discussion of the relevance of these results for state of the art 
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experiments. Recently, multiple scattering experiments [6, 26] were performed on disordered 
media of titania powders with thicknesses from 5 to 20 ;um, and a mean free path d, ^ 0.9 
/im. For an average refractive index n ^ 1.34 [26], this implies Thouless times tu = (c/n)£/3 
ranging from 3.8 x 10~^^s to 6.0 x lO^^^s, with c the speed of light in vacuum. Thus, photons 
with subpicosecond or subnanosecond coherence time will allow us to explore both scenarios 
outlined above, t^ ^ to and t^ ^ Id-, and to quantify photon entanglement as well as the 
transport characteristics of the medium. 
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